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Abstract 

We study the rigidity of body- and- cad frameworks which capture the majority of the geometric con- 
straints used in 3D mechanical engineering CAD software. We present a combinatorial characterization 
of the generic minimal rigidity of a subset of body-and-cad frameworks in which we treat 20 of the 21 
body-and-cad constraints, omitting only point-point coincidences. While the handful of classical combi- 
natorial characterizations of rigidity focus on distance constraints between points, this is the first result 
simultaneously addressing coincidence, angular, and distance constraints. Our result is stated in terms 
of the partitioning of a graph into edge-disjoint spanning trees. This combinatorial approach provides 
the theoretical basis for the development of deterministic algorithms (that will not depend on numerical 
methods) for analyzing the rigidity of body-and-cad frameworks. 



1 Introduction 

We study ftody-and-cad frameworks composed of rigid bodies with pairwise coincidence, angular and distance 
constraints placed between geometric elements (points, lines, planes) rigidly affixed to each body. A frame- 
work is flexible if the bodies can move relative to each other while respecting the constraints; otherwise. 



it is rigid (see Figures la and lb). A fundamental problem for both the user and the CAD software is to 



determine whether a system is flexible or rigid. 

Contributions. In this paper, we present a combinatorial characterization of the generic minimal rigidity 
of a subset of body-and-cad frameworks: we treat 20 of the 21 constraints, omitting point-point coincidences. 
The equations governing point-point coincidences exhibit special algebraic behavior (see Appendix [b|) . 

Body-and-cad frameworks on n bodies without point-point coincidences can be modeled by a graph called 
a (fc, (/)-frame. Our main theorem states that generic minimal rigidity of a (fc,(7)-frame is equivalent to a 
partitioning of the graph's edge set such that each partition consists of the edge-disjoint union of 3 spanning 
trees. Combinatorial characterizations of rigidity typically lead to quadratic time algorithms, and we expect 
this result to have similar implications for the development of efficient algorithms. 

Significance. Among the many classes of structures whose rigidity is studied, there are only a handful of 
cases for which combinatorial characterizations have been proven, with Laman's 2D bar-and-joint [TOj and 
Tay's body-and-bar [20. results the most well-known. Furthermore, classical work has focused mainly on 
distance constraints. To the best of our knowledge, ours is the first combinatorial rigidity characterization 
that simultaneously treats angular and other constraints. 

Motivation. In the 3D "assembly" environment of the popular CAD software SolidWorks users place 
"mates" among parts (rigid bodies). The mates specify constraints (e.g., distance, angular, coincidences) 



* Final version to appear in Symposium on Solid and Physical Modeling '12 and associated special issue of Computer Aided 
Design. 
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(a) A minimally rigid body-and-cad framework with (i) a 
line-line coincidence along the dotted blue axis, (ii) a point- 
plane coincidence between the red point r and shaded plane, 
and (iii) a point-point distance between the green points p 
and q. 



(b) Removing the point-point distance constraint (iii) re- 
sults in a flexible framework with 1 degree of freedom (ro- 
tation about the shared axis). 



Figure 1: A 3D body-and-cad framework. 



among geometric elements (e.g., points, lines, planes) identified on the parts. Informative feedback is im- 
portant as explicit visualization of mates is difficult, and mechanical engineers are usually not experts in 
reasoning about the underlying geometric constraint systems. 

One of the most difficult situations arises when a new mate is specified that causes an inconsistency in 
the design. As depicted in Figure [2] SolidWorks gives feedback in the form of an alert window and highlights 
a list of previous mates "overdefining" the system. This list is intended to assist the user in resolving 
the inconsistency, but is often too large (sometimes, the entire set of mates) to be helprooful. Current 
software cannot reliably detect a minimally dependent set of constraints. Indeed, the system in Figure [2] 
reveals the severity of software limitations, as the lift is essentially 2-dimensional. It is composed of two 
"scissor" parts attached by struts connecting pairs of matching joints. Since the "scissor" parts must move 
identically, analysis can be reduced to considering a single "scissor," which itself is a planar mechanism. The 
development of a complete rigidity theory would provide a rigorous mathematical understanding of these 
geometric systems. 
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Figure 2: When a new mate is added that results in an inconsistency, previous mates "over defining the 
assembly" are highlighted. However, the highlighted set is not necessarily minimal. SolidWorks model from 
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Potential for algorithmic impact. Given the pebble game algorithms for sparsity in 'TTj, we anticipate 
our result will lead to efficient approaches for determining the rigidity of body-and-cad frameworks as well 
as finding rigid components and detecting minimally dependent sets of constraints. Algorithms based on our 
characterization would provide the capability of giving valuable feedback to users designing highly complicated 
CAD systems. 

Related work. Geometric constraint systems are at the core of constraint-based CAD software and have 
been studied from several perspectives. 

In the CAD community, emphasis has been on finding solutions for realizing geometric constraint sys- 
tems. A common approach is to use a decomposition-recombination scheme based on graph algorithms and 
numerical solvers; see, e.g., [9 19 for a survey of standard techniques. These decomposition schemes often 



rely on breaking a system into rigid (also called "well-constrained") sub-systems. Difficulties arise even in 
studying small rigid sub-systems. The work of Gao et al. [4] enumerates and analyzes all "basic configura- 
tions" with up to six geometric primitives, then presents a method for finding their solutions; the authors 
observe that the most challenging configurations included constraints involving lines (in contrast to those 
involving only points or planes). 

In classical rigidity theory, bar-and-joint structures composed of universal joints connected by fixed length 
bars, i.e., point-point distance constraints, are studied. For an overview of combinatorial rigidity theory, see 
the texts |6, 7 . 

Combinatorial properties of rigidity are usually tied to sparsity counting conditions. A graph G is {k,i)- 
sparse if the induced subgraph on any subset of n' vertices contains at most kn' — £ edges. Under certain 
conditions, there are efficient algorithms to detect (fc, £)-sparsity. If G is (fc, ^)-sparse and < £ < 2k, then 
its edges give rise to a matroid. Within this matroidal range of values of k and £, (fc, £)-sparsity is determined 
by a family of pebble game algorithms which run in 0{n^) worst case time 
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In 2D, bar-and-joint rigidity is characterized by Laman's theorem with (2, 3)-sparsity; other characteriza- 
tions followed, some based on decompositions of the edge set into trees [2| [T3|[l5] . However, 3D bar-and-joint 
rigidity is not well-understood, and a combinatorial characterization is arguably the biggest open problem in 
rigidity theory. While (3, 6)-sparsity is necessary, it is not sufficient. Figure |3] depicts the classical "double 
banana" counterexample; although it satisfies the necessary counts, the structure is flexible, as the "bananas" 
can rotate about the dotted axis. 




Figure 3: Double banana counterexample shows that (3, 6)-sparsity is not sufficient for 3-dimensional bar- 
and-joint rigidity. The structure is flexible, as the "bananas" joined at points a and b can rotate about the 
dotted axis. 



While bar-and-joint structures in 3D are not well-understood, a closely related structure called the body- 
and-bar structure is. The body-and-cad constraints that we consider include point-point distance constraints, 
which are simply bars - the body-and-bar rigidity model is a special case of body-and-cad. Body-and-bar 



rigidity was characterized by Tay 20 ; an alternate proof was given by White and Whiteley j22j. We 
generalize White and Whiteley's proof to obtain our result. Tay's theorem [20] states that d-dimensional 
body-and-bar rigidity is equivalent to an associated graph consisting of ('^^^) edge-disjoint spanning trees; 
as a consequence of Nash- Williams and Tutte's theorems [l4{[2T] on the equivalence of edge-disjoint spanning 
trees and sparsity, d-dimensional body-and-bar rigidity is characterized by ((''2^), (''2^))-sparsity. 

Other constraints motivated by CAD applications have been studied in rigidity theory. Servatius and 
Whiteley give a combinatorial characterization of rigidity for frameworks with direction (orientation of 
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the vector defined by two points with respect to a global coordinate frame) and distance constraints 18 
However, angular constraints among points, even in the plane, have proven very challenging; in fl6l Saliola 
and Whiteley showed that deciding independence of circle intersection angles in the plane has the same 
complexity as determining 3D bar-and-joint rigidity. Angular constraints between lines and planes or rigid 



bodies, though, is well-understood 12 



The results from the rigidity theory community have not been directly applied in CAD research partly 
because rigidity for CAD systems has not yet been fully characterized, and partly because of issues related 
to genericity (discussed in detail in Section 5.11. Note that if a framework is rigid, then the set of all rigid 
frameworks with the same combinatorics is an open subset of its parameter space. Michelucci and Foufou [I] 
use this observation and present probabilistic algorithms which rely on analyzing a generic "witness" which 
shares the incidence constraints of the original system, but may have different parameters for the distance 
and angular constraints. Even in the plane, where Laman's theorem characterizes bar-and-joint rigidity, the 
restriction of genericity poses a subtle challenge. The work of Jermann et al 23 observes that coincidence 
constraints, when expressed as distance constraints with value 0, exhibit non-generic behavior; a similar 
situation arises for parallel constraints. This is addressed through a notion of "extended structural rigidity" ; 
Jermann et al make assumptions similar to the genericity conditions for the characterization we present in 
this paper. 

The foundations of infinitesimal body- and- cad rigidity theory were introduced in From the develop- 
ment of a rigidity matrix, a natural nested sparsity condition was identified as a necessary, but not sufficient 
condition for rigidity. The example given there highlighted the need for a better understanding of angular 
constraints. The need for special treatment of angular constraints was implicitly observed in [s], but, to the 
best of our knowledge, there has been no other work providing an explicit understanding of their behav- 
ior. With both the flavor of 3-dimensional bar-and-joint and body-and-bar rigidity, the question remained 
as to where the difficulty of body-and-cad rigidity lay. The results we present here provide an answer for 
the majority of body-and-cad constraints, but indicate that point-point coincidences may pose a significant 
challenge (see Section [5^ . 

Structure. In Section [2| we provide the preliminary concepts necessary for understanding the foundations 
of infinitesimal body-and-cad rigidity and give an example to highlight the distinct behavior of angular 
constraints. In Section [3j we develop the rigidity matrix for a more general combinatorial object called a 
{k, (7)-frame, which is used to prove our main theorem in Section|4] In Sectionjsj we discuss the subtleties and 
limitations resulting from genericity assumptions as well as the challenges posed by point-point coincidence 
constraints. We conclude with future work in Section [6l 



2 Preliminaries 

For a geometric constraint model, there are three levels of rigidity theory: algebraic, infinitesimal and 
combinatorial. In algebraic rigidity theory, a system of equations expressing the geometric constraints 
is studied; a solution to this system corresponds to a realization of the structurej^ In infinitesimal rigidity 
theory, the first-order behavior of the (usually quadratic) system of equations expresses the constraints in 
terms of instantaneous motions and a rigidity matrix. In combinatorial rigidity theory, rigidity is defined 
terms of the rank of this rigidity matrix, so that generically rigidity depends only on properties of a certain 
graph. In this paper, we are concerned only with infinitesimal and combinatorial body-and-cad rigidity. 

The infinitesimal body-and-cad rigidity theory was first presented in |8 and relies on expressing 
constraints in the Grassmann-Cayley algebra through the rows of a rigidity matrix. For completeness, we 
provide an overview of the foundations developed in [sj. 

Body-and-cad constraints. The body-and-cad frameworks that arise from CAD software are naturally 
3-dimensional. There are 21 coincidence, angular and distance constraints that can be placed between a pair 
of geometric elements (points, lines or planes) on rigid bodies, which we enumerate below for clarity: 

• Point-point constraints. Coincidence, distance. 

^In Appendix [c] wc provide the formal definitions associated with body-and-cad algebraic rigidity theory, including the 
family of "length" functions required to describe the geometry of the system. 
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• Point-line constraints. Coincidence, distance. 

• Point-plane constraints. Coincidence, distance. 

• Line-line constraints. Parallel, perpendicular, fixed angle, coincidence, distance. 

• Line-plane constraints. Parallel, perpendicular, fixed angle, coincidence, distance. 

• Plane-plane constraints. Parallel, perpendicular, fixed angle, coincidence, distance. 

Infinitesimal constraint equations. To derive equations from the geometric constraints, we must first 
consider instantaneous rigid body motion in M^. As a consequence of Chasles' Theorem (see, e.g., 17 ), any 



instantaneous rigid body motion may be described by a twist (translation and rotation about a specified 
twist axis), which itself is represented by a 6- vector s = (a;, v). The 3- vector a; describes the angular velocity: 
the direction of the twist axis and rotational speed about it. The 3- vector v can be used to decode the rest 
of the twist axis and translational speed along it. 

A primitive constraint between two bodies i and j is encoded by a single homogeneous linear equation 
on the twists s^ = {uJi,Vi) and Sj = {u}j,Vj). Each body-and-cad constraint is associated to a number 
of primitive constraints (which intuitively affect at most one degree of freedomjj^ A distinction is made 
between primitive angular and blind constraints: a primitive angular constraint may affect only a rotational 
degree of freedom, while a primitive blind constraint may affect either a rotational or translational degree of 
freedom. Equations corresponding to angular constraints have zeroes in the coordinates corresponding to v. 

The main contribution of [8j was the algebro-geometric derivation of these equations which are collected 
together in the rigidity matrix to express the infinitesimal body-and-cad constraints. Given a structure with 
n bodies, the matrix has 6n columns; these are arranged so that the 6 columns for the ith vertex correspond 
to (vi, — Wijj^ For the combinatorial characterization that we present in this paper, we are only 
concerned with the pattern of non-zero entries in the rigidity matrix. 

Elements of the kernel of the rigidity matrix may be interpreted as the infinitesimal motions of the 
framework. If the only infinitesimal motions are trivial, i.e., assign the same twist to each body, the framework 
is infinitesimally rigid; otherwise, it is infinitesimally flexible. Since we will only work in the infinitesimal 
rigidity theory, for brevity, we will drop "infinitesimally" for the remainder of this paper. A minimally rigid 
framework is one that is rigid, but becomes flexible after the removal of any (primitive) constraint. 

Example. We provide a small example to illustrate the concepts of primitive constraints and their further 
separation into angular and blind constraints. To emphasize that we are only concerned with the pattern 
of zero and non-zero entries, we use *-entries to indicate values that are generically non-zero. 



The body-and-cad framework depicted in Figure la is composed of two bodies (a cube and hexagonal 



prism) sharing three constraints: (i) a line-line coincidence, (ii) a point-plane coincidence and (iii) a point- 
point distance. This example is minimally rigid; for instance, removal of the (primitive blind) point-point 



distance constraint results in a flexible structure with one degree of freedom (see Figure lb ) . 

The constraints for this framework are inflnitesimally expressed via a rigidity matrix with the pattern 
depicted in Figure |4j Observe that the line-line coincidence constraint is associated to four rows of the 
matrix; this corresponds to four primitive constraints. In fact, rows 1 and 2 express a line-line parallel 
constraint and have values in three columns associated with each body; each of these rows is expressing a 
primitive angular constraint. The remaining four rows express primitive blind constraints. 

Novelty. We emphasize the two main differences between the development of the infinitesimal rigidity 
theory for the body-and-cad (S) and body-and-bar 20 22] models. In the rigidity matrix associated to a 



body-and-bar framework, each row corresponds to a bar, or equivalently, a distance constraint between 
two bodies. By contrast, in order to encode a single constraint in a body-and-cad framework, we may 
need multiple rows in the rigidity matrix; the concept of primitive constraints addresses this by associating 
exactly one constraint with each row. 

The second difference appears in the separate treatment of primitive angular and blind constraints, 
represented by red and black edges, respectively, in a primitive cad graph, defined below. This was the 



"^Four basic constraints, described in more detail in Appendix|A] expressed 20 of the constraints; in Appendix|B] we discuss 
why we cannot address point-point constraints with our proof technique. 

^The re-ordering and negation are technicaUties arising from the development of the rigidity matrix. 



5 



Vl 



-(JJl 



V2 



-OJ2 



line- line 




* * * 


000 


* * * ] 


coincidence 


* ^ * 


* * * 


* * * 


* * * 




* * * 


* * * 


* * * 


* * % 


point-plane distance 


* * * 


* * * 


* * * 


* * * 


point-point distance 


* * % 


* * * 


* * * 


* * * 



Figure 4: The pattern of the rigidity matrix for the structure depicted in Figure la 



main challenge left open by the natural nested sparsity condition identified as a necessary, but not sufficient, 
condition for body-and-cad rigidity. The example given on page 27 of [s] highlighted the need for a better 
understanding of the interplay of angular and blind constraints. 



2.1 The primitive cad graph 

The infinitesimal theory of ^ lets us associate a primitive cad graph H — [V^BU R) to each body-and-cad 
framework. Here, H is a multigraph with a vertex for each body and an edge for each primitive constraint. 
Moreover, the edges are partitioned into two sets, a set B of black edges and a set R of red edges, representing 
primitive blind and angular constraints, respectively. See Figure [5] for the primitive cad graph associated to 
the minimally rigid framework in Figure [ij 



line-lini 
coincidenci 




pom'Mjomt 
coincidence distance 




Figure 5: The primitive cad graph associated to the framework from Figure la has 2 red (dashed) edges, 
representing primitive angular constraints, and 4 black (solid) edges, representing primitive blind constraints. 

Recall that each primitive constraint is expressed as a single linear equation. In fact, this equation can 
be encoded by a 6- vector whose entries are derived from the specific geometry of the framework. We define a 
primitive cad framework H{p) to be the graph H in which each edge is labeled with this 6-vector. Let m be 
the number of edges in H and mji — \R\ be the number of red edges. Our main theorem for 3-dimensional 
body-and-cad rigidity is stated in terms of the primitive cad graph H. 

Theorem 1. A 3- dimensional body-and-cad framework (in which no point-point coincidence constraints are 
present) is generically minimally rigid if and only if in its associated primitive cad graph H = {V, RUB), 
there is some set of black edges B' Q B such that 

1. B\B' is the edge-disjoint union of 3 spanning trees, and 

2. RU B' is the edge-disjoint union of 3 spanning trees. 

Body-and-cad frameworks are a special case of the more general (fc, g)-bsimes introduced in Section |3] 
Theorem [1] follows from a result (Theorem 16 1 in this more general setting, which we prove in Section |4] 



3 From body-and-cad frameworks to {k, g)-fraines 



We begin in Section [3.1| by introducing the (/c, g)-frame, which generalizes both the notion of a body-and-cad 
framework and the fc-frame of |22|, and its associated rigidity matrix. We then present a pure condition 
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for the minimal rigidity of a (fc, g)-frame in Section 3.2 From the point of view of apphcations, our resuh 
is most interesting for 3-dimensional body-and-cad frameworks; however, the techniques do not depend on 
dimension, so we work in the most general setting. 



3.1 (fc, (7)-frames 

The proof of our main result depends only on the pattern of zeroes in the rigidity matrix of a body-and-cad 
framework, and not on their derivation or interpretation. Thus, we may view the pattern in the rigidity 
matrix associated to a primitive cad framework as an instance of a rigidity matrix associated to a more 
general (/c, (7)-frame, which we introduce below. We begin with some preliminary combinatorial definitions. 

Definition 2. A multigraph H ~ {V, E = B U R) is a bi-colored graph, where V ^ {I, . . . ,n} is the vertex 
set, and E is a set of m edges, nig of which are black andmn of which are red. For brevity, we use "graphs" 
to refer to multigraphs. 

We generalize the notion of a fc- frame in Definition 1.7 of [22] . 

Definition 3. A (fc,(7)-frame H{p) is a bi-colored graph H = (V, E = B U R) together with a function 
p : E ^ M*^, where p{r)j = for j = 1, . . . , k ~ g and r E R. 

Then, primitive cad frames in 3D are simply (6, 3)-frames. The /c-frames of White and Whiteley are 
(k, g)-bames in which R is empty (hence the value oi g is irrelevant). 

Associated to each (fc, 5)-frame is a rigidity matrix which we define below. 

Definition 4. Given a (fc, g)-frame H[p), we define a matrix M{H{p)) with k columns per vertex. For each 
edge e = uv with u < v, define a row with p(e) in the fc columns for u and — p(e) in the fc columns for v. 

Motions assign a vector of length fc to each body; for structures in dimension d, k ^ (''+^), and each 
vector represents a compatible infinitesimal motion. 

Definition 5. A motion of a (k,g) -frame is a vector of length kn that is orthogonal to the row space of 
AI{H{p)). A trivial motion assigns the same k-vector to each body, i.e., has n copies of the same k-vector. 

We are concerned with minimal rigidity, i.e., structures with no redundant constraints. 

Definition 6. A {k, g)-frame H{p) is minimally rigicj^ if it becomes flexible after the removal of any edge 
from H . 

Remark 7. _De/iTOtzons and sei the pattern of zeroes in the matrix M{H{p)). However, the rows of 
the body-and-cad rigidity matrix presented in |8j are more specialized. In Appendix^^we give an alternative 
derivation of equations for two of the basic body-and-cad constraints so that the equations conform to the 
pattern we have set. 



3.2 The pure condition for minimal rigidity 

As with other classical rigidity results, our characterization relies solely on the combinatorics of a body-and- 
cad framework, which applies generically. In this section we introduce the pure condition for the minimal 
rigidity of a (fc, (7)-frame, generalizing the setup in 22 . In fact, the existence of a pure condition is what 



justifies our genericity assumptions, which we discuss in more detail in Section 5.1 



The definition below generalizes necessary edge-counting conditions for minimal rigidity, taking angular 
constraints into account. 

Definition 8. Let fc and g be positive integers. A bi-colored graph H is (k, g)-comited if 
1. m — kn — k (count on total edges), and 
^- ™fl ^ giT' ^ g (count on red edges). 



Also called "isostatic" in the literature. 
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The matrix M{H{p)) of a (fc, f;)-counted {k, g)-bame is not square (it has kn~ k rows and kn columns). 
It is straightforward to observe that any minimally rigid (k, g)-bame must have kn—k rows since the trivial 
motions will always be in the kernel of M{H{p)). If we fix or "tie down" one of the bodies in the framework, 
then these trivial motions are no longer allowed. 

Definition 9. The basic tie-down of a (k, g)-counted {k,g)-frame is a k x kn matrix T{k) with the identity 
matrix in the first k columns and Os everywhere else. 

Given a [k, g)-counted (k, g) -frame H{p), the (square) matrix Mt{H{p)) isAI{H{p)) withT(k) appended 
to the bottom. 

Thus, the rigidity of a (fc, g)-frame H{p) satisfying the necessary counts may be determined by computing 
the determinant of the matrix Mt{H{p)). 

Proposition 10 (Proposition 2.5 from |22j). A [k, g)-counted {k, g)-frame H{p) is minimally rigid if and 
onlytfMT{H{p))^0. 

In fact, since we would like to understand rigidity more generally than for a single set of values for p, 
we would like to think of Mt{H{x.)) as a polynomial in indeterminates. As in 22], we make the following 
definition. 

Definition 11. Given a (k,g)-frame H{p), we define the generic {k, g)-frame H{x.) on the underlying labeled 
graph H by setting all generically non-zero entries of the k-vectors labeling the edges to be algebraically 
independent indeterminates x(6)j for j — I, . . . ,k, and x(r)j for j = fc — g + 1, . . . , fc, where b d B and 
r e R. 

We define the pure condition to be Ch{^) = det Mt{H{x)) so that C/f(x) is a polynomial in iV = 
krriB + gmn variables. Generic rigidity is expressed in terms of the pure condition. 

Definition 12. A [k, g)-counted {k,g)-frame i?(x) is generically minimally rigid if there exists a function 
p : — > M'^ for which Ch{v) 7^ satisfying p{r)i ~ for i — 1, . . . ,k — g for all r £ R. 

We require one final combinatorial concept for analyzing the pure condition. 

Definition 13. A {k,g)-ia,Ti (j) of a [k,g)-counted bicolored graph is a partitioning of E into n — 1 ordered 
sets {(p2, ■ • . , (pn) such that 

1. each (pi — (0i_i, . . . , 0i.fe) contains exactly k edges incident to vertex i, and 

2. each <f>i contains < g (red) edges from R. 

We denote the determinant of the matrix whose rows are x((/)j^i), . . . ,ii.((f)i_k) by [pi]. 

Definition 14. Two {k, g)-fans p and p' are called distinct if there exists a vertex i such that pi ^ p'i, as 
unordered sets. 

We can represent distinct (fc,(7)-fans by orienting H: if edge e ^ pi, orient it so that its tail is at vertex 
i. Then, given a {k,g)-ia,n {p2, . . . ,pn), its {k,g)-fan diagram is the oriented multigraph F = (y,A), where 
ij € ^ if and only if ij e pi. Note that the same {k,g)-ia,n diagram will represent multiple (k, g)-fans, as 



depicted in Figure 6b but exactly one distinct (k, g)-fan. 

As sl consequence of Proposition 2.12 from [22], we obtain the following. 

Proposition 15. Ch{x) = ± [^2] • • • [pn], for all distinct {k,g)-counted {k,g)-fans p of H . 

Proof. If we compute det AIt{H{x.)) via a Laplace expansion on the first k columns of Mt{H{x.)), we see 
that the determinant is actually the determinant of the {kn — fc) x {kn — k) matrix formed by deleting the 
first k columns and last k rows. We can in turn compute the determinant of this submatrix via a Laplace 
expansion taking the k columns of each vertex at a time. This requires that we compute the determinant of 
matrices formed by partitioning the rows of our matrix so that k rows are associated to each vertex 2, . . . , n. 
Each row corresponds to an edge, and if we have a submatrix containing more than g rows associated to 
red edges, then the corresponding determinant is zero. Therefore, C/f (x) is the sum of terms of the form 
±[(/)2] • • • [pn], where p — {p2, ■ ■ ■ , Pn) is a distinct (fc, g)-fan. The sign of each term is determined by the 
rules for Laplace expansion, the permutation p, and whether the vectors used in the Laplace expansion are 
x(e) or — x(e). □ 
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(a) Example (3, l)-countod bicolored 
graph has 9 total edges, 2 of which are 
red. 




(b) This (3, l)-fan diagram corre- 
sponds to multiple (3, l)-fans, in- 
cluding: ({c,d,e),{g,h,i),(a,b,f)) and 
((ex,d),(i,g,h),(f,a,b)). 




(c) Black edges without bolded edge 
form 2 edge-disjoint spanning trees 
(dotted and dashed); red edges with 
bolded edge form a single spanning tree 
(solid). 



Figure 6: An example (3, l)-couiited graph satisfying the conditions of Theorem 
edges: c between vertices 1 and 2, and h between vertices 3 and 4. 
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The graph has two red 



4 A combinatorial characterization of minimal rigidity 

Given the setup of ^ and the previous sections, our main resuh may be stated and proved in a purely 



combinatorial setting. Theorem 16 and its proof are a generalization of the body-and-bar characterization 
theorem of White and Whiteley. 

Theorem 16. A (k, g)-counted (k, g)-frame -ff(x) is generically minimally rigid if and only if there exists 
some set of black edges B' C B such that 

1. B\B' is the edge-disjoint union of k ^ g spanning trees, and 

2. RU B' is the edge-disjoint union of g spanning trees. 



Figure 6c depicts an example of a (3, l)-counted bi-colored graph satisfying the conditions of Theorem |16[ 



Proof. (=>) Assume Ch{x) 0. Consider a Laplace expansion along the last k rows of Mt{H{'x.)). Since 
the only non-zero entries in these rows occur in the k x k identity matrix appearing in the tie-down in the 
first k columns, detMy(i?(x)) — detA, where A is the submatrix of Mt{H(x.)) formed by the first kn — k 
rows and last kn — k columns. Since C//(x) — det Mt{H{x.)), then detA ^ 0. 

Now consider a Laplace expansion of det A using (n — 1) x (n — l)-minors so that det A is a sum of terms 
II = det Ai ■ ■ ■ det Ak, where for j = 1, . . . , fc, is an (n — 1) x (n — 1) submatrix of A using the jth columns 
associated to each of the n — 1 remaining vertices and some choice of n — 1 rows. If det A is non-zero, then 
some term fi = det Ai ■ ■ ■ det A^. is non-zero, which implies that det Ai, . . . , det A^ are all non-zero. Each 
submatrix Aj has one column per vertex, and the rows of Aj are just the rows of the incidence matrix of an 
orientation on the edges of H multiplied by non-zero scalars. Since det Aj is non-zero, Aj is the incidence 
matrix of a subgraph of H with n — 1 edges and no cycles. Hence, Aj describes a spanning tree Tj of H. 

If r is a red edge, then x(r)j = for j ^ 1, . . . , k — g. Since det Aj is non-zero for all j, the edge r must 
be in one of the trees T^-g+i,. . . ,Tj.. Let B' C B he the set of black edges in T^^g+i, . . . , T^. Then B\ B' 
is the edge-disjoint union of Ti, . . . , Tfc-g and RU B' is the edge-disjoint union of the g remaining spanning 
trees. 

(■^=) Let -ff (x) be the generic (fc, g)-ira,me, where H = {V, E = B Li R) is a, {k, gj-counted graph, and let 
S' C _B be a set of edges such that 

1. B \ B' is the edge-disjoint union of fc — g spanning trees Ti, . . . , Tk-g, 

2. RUB' is the edge-disjoint union of g spanning trees Tfc_g+i, . . . , T^. 
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We show that dct Mt{H{-x)) is not identically zero. To do this it suffices to show that there exists a function 



x.' : E M.'' such that det AIt(H(x')) 7^ 0. For j — 1, . . . , k and i — 1, . . . , k, let aj i denote an indeterminate. 



For each j and e G Tj let 

x'(e) 



(aja,a-,,2, . . • ,aj,fc) if e G B \ B' 

(0,0, . . . ,0,aj,fe_g+i, . . . ,aj,fc) if e e i?U 



By Proposition 15 Ch(x) — ± n"^2[0i]- We claim that the only non-zero term in the pure condition 
stems from a single distinct (/c, g)-fan (j). Root each tree at vertex 1. Each tree has a unique path from vertex 
i to the root, so we can orient all edges toward the root. Define (t>i to be the ordered set of all edges incident 
to vertex i pointing toward the root, where the j'th edge in (pi is the edge in Tj. By construction, each (pi 
has at most g red edges, so [cpi] ^ 0. 

Consider another distinct (fc,(7)-fan 0' ^ 0. There must be some vertex i such that contains at least 
two edges ei and 62 from the same tree Tj. Hence, x'(ei) = x'(e2) and [i/)^] — 0, and we conclude that </> is 
the only fan for which n"^2['/'i] is non-zero. □ 

We now prove Theorem [l] which gives a combinatorial characterization for 3-dimensional body-and-cad 
rigidity (omitting point-point coincidence constraints). 

Proof of Theorem^ The result follows from Theorem [16] as a primitive cad graph for a body-and-cad frame- 
work with no point-point coincidence constraints is a (fc, g)-framc with fc = 6 and .g = 3. □ 



5 Considerations and challenges 

Before concluding, wc discuss practical implications of the genericity assumptions that result from the pure 
condition, then provide context as to why point-point coincidence constraints pose a challenge. 

5.1 Genericity 

Like all combinatorial characterizations of rigidity, our main theorem gives a condition for generic minimal 
rigidity. This notion of genericity is subtle as there are currently no geometric conditions implying or 
characterizing it (even in the most well-studied bar-and-joint rigidity model). Common assumptions in 
computational geometry, such as the general position of joint coordinates, are not strong enough to ensure 
the genericity of a bar-and-joint structure, and we have additional types of constraints to consider. Thus it 
is difficult to test when a given framework satisfies the genericity assumptions needed for the combinatorial 
analysis. 



We defined generic rigidity using the pure condition (see Definition 12 ) and now discuss why this definition 
is appropriate. Recall that = kms + grriR. The set of all points q € with Ci/(q) = is a closed 
subset that either has dimension less than N or is all of R''^. Consequently, the set of points q € such 
that Cjy(q) ^ is open. This implies that, if there exists one q G such that Ch{<i) 7^ 0, then the pure 
condition is non-zero for all points in some small neighborhood of q. Moreover, if C^f(q) ^ for some q, 
then the set of points with Cjy(q) = is a proper closed subset of with measure zero, and C/f (q') 7^ 
generically, i.e., for almost all values q' G M.^ . 

Therefore, Definition [12] can be viewed as an extension of Tay's notion of generic rigidity |20j to {k,g)- 
frames. A (A;, (7)-frame -ff(x) is generically minimally rigid if and only if C/f(x) is not identically zero. If, 
however, Cij(x) is not the zero polynomial and Ch{p) = for some p G M^, then H{p) is (infinitesimally) 
flexible. We refer to p as a non-generic point, and say that we have a non-generic embedding of the {k,g)- 
frame H{p). If the polynomial C/f (x) is identically zero, then i?(x) is generically flexible. 

A non-generic framework: Pappus's Theorem 

As observed by Jermann et al ,23 , there are inherent problems in analyzing the geometry of a specific 
realization of a framework using only the combinatorics of a graph, as a particular embedding or framework 
may fail to be generic. 

This problem is straightforward to observe in the analysis of body-and-bar frameworks, in which all 
constraints are point-point distance constraints. For example, consider a graph with vertices vi and V2 
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joined by 6 edges; this is the combinatorial model for two rigid bodies joined by 6 bars. If the bars are 
attached at points that are chosen generically, then the framework is rigid, matching the combinatorial 
analysis provided by Tay's Theorem |20j. However, if the points are not sufficiently generic - for example if 






(a) The 6 points A, B, C, a, b, c in the non-generic position (b) Circular vertices are associated to the bodies from points 
from the hypothesis of Pappus's Theorem. and square vertices are associated to the bodies from lines. 

Vertices are connected if there is a point- line coincidence be- 
tween the corresponding rigid bodies. There are two edges 
for each point-line constraint in the primitive cad graph, 
but we have drawn only one to aid readability. The vertex 
labeled n line represents the line PabPbc- 

Figure 7: The point-line coincidences described by Pappus's Theorem. 

all 6 bars attach to a single point on one of the bodies - then the framework is flexible; for this non-generic 
embedding, Tay's theorem does not apply. Intuitively, "reusing" the same point results in a dependency 
that is not encountered in the generic case. 

The constraints in a body-and-cad framework present similar challenges in applying combinatorial analy- 
sis. To illustrate this, we present a body-and-cad framework along with a non-generic embedding, modeling 
Pappus' Theorem, in which a generically independent constraint becomes dependent. 

Theorem 17 (Pappus). (See Figu re \ 7a[ ) Let A,B,C and a,b,c be two sets of 3 collinear points in the 
plane. Then the 3 points pab = Ah H Ba, pac — Ac H Ca, and pbc — Be H Ch are collinear. 

To model Pappus's Theorem, we use a framework with 18 bodies: 9 to model the points, and another 9 
to model the two lines ABC and a6c, the six lines Ah, Ac, Ba, Be, Ca, Ch, and the line PatPbc- If we encode 
all of the point-line coincidences depicted in diagram (a) in Figure [t] then the point-line coincidence between 
Pac and PabPbc IS implied by the other constraints as a consequence of Pappus's Theorem. Thus, there is 
a dependency among the geometric constraints which should correspond to a dependency in the rigidity 
matrix. 

However, a combinatorial analysis of the associated primitive cad graph (depicted schematically in Figure 
7b) based on Theorem [I] indicates that the constraints are generically independent. Indeed, we can create a 
generic embedding of the same primitive cad graph in which each constraint uses unique geometric elements: 
for each point-line coincidence, identify a distinct point and distinct line on the pair of constrained bodies. 
If the 27 points and 27 lines are sufficiently generic, then the associated rigidity matrix has full rank, and 
the constraints are generically independent. 



5.2 Point-point coincidence constraints 

The result given in this paper does not address point-point coincidence constraints. In addition to the alge- 
braic difficulties (described in Appendix [b|) that prohibit the inclusion of point-point coincidence constraints 
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in the definition of a (fc, (7)-franie, point-point coincidences seem inherently more challenging to address. 

It may be that characterizing 3D point-point coincidences is just as challenging as understanding 3D 
bar-and-joint rigidity. Indeed, we can construct an analogue of the classical "double banana" example with 
a body-and-cad framework that has two rigid bodies joined by two point-point coincidence constraints. The 
two rigid bodies are free to rotate about the line joining their two points of intersection (labeled a and b in 
Figure |3|, so the body-and-cad framework is flexible. 

In this example, each point-point coincidence constraint corresponds to three primitive blind constraints 
(see Appendix|B]) , and the primitive cad graph contains 2 vertices and 6 black edges between them. A natural 
extension of Theorem [T] would characterize generic minimal rigidity with a partitioning of the edge set into 
6 edge-disjoint spanning trees (there are no red edges). While this property is necessary, it is not sufficient. 
The primitive cad graph for the "double banana" satisfies the condition, but is fiexible, thus serving as a 
counterexample. Therefore, a combinatorial characterization of body-and-cad frameworks including point- 
point coincidences must differ from the characterization for body-and-cad frameworks presented here. 



6 Conclusions and Future Work 

We have presented a combinatorial characterization of body-and-cad rigidity of structures using 20 of the 
21 pairwise constraints; we exclude point-point coincidence constraints. This characterization is expressed 
in terms of the edge set of an associated graph being partitioned into two sets of 3 edge-disjoint spanning 
trees. 



Future work and open problems. The discussion in Section 5.1 highlights a practical, but notoriously 
difficult, problem in rigidity theory: determine geometric necessary or sufficient conditions for a specific 
embedding of a body-and-cad framework to be generic. 

We suspect that analyzing point-point coincidences in body-and-cad frameworks may be just as hard as 
analyzing 3D bar-and-joint rigidity. Understanding precisely where the challenge lies merits a comprehensive 
investigation, and the relationship between 2D body-and-cad and 2D bar-and-joint structures is a natural 
starting point. 

Algorithms. The combinatorial property of Theorem [T6| must be matroidal as it characterizes the indepen- 
dence of rows in a matrix. Due to the intimate relation between edge-disjoint spanning trees and sparsity 
counts, we expect to obtain efficient algorithms for body-and-cad rigidity by generalizing the pebble game 
algorithms. Such algorithms would not only decide whether a given framework is rigid or flexible, but 
should be able to detect rigid components and circuits, minimally dependent sets of (primitive) con- 
straints. This would provide valuable feedback to CAD users when the addition of a constraint causes an 
inconsistency. 
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A Alternative development of body-and-cad rigidity matrix 



In [s], a body-and-cad rigidity matrix was developed by expressing each of 20 possible 3-diniensional con- 
straints in terms of four basic constraints: 
(i) basic line-line non-parallel fixed angular, 
(a) basic line-line parallel, 
(Hi) basic blind orthogonality, and 
(iv) basic blind parallel. 

Point-point coincidence constraints required a separate development, and the resulting rows do not have the 
pattern of generically non-zero entries required for the (fc, g)-baines in this paper. 

Constraints (i) and (ii) are angular constraints, as the corresponding twists have v = 0. Each basic 
constraint corresponds to either one or two rows in the rigidity matrix. 

The basic angular constraint (i) requires that the angle between two vectors a and b be maintained. 
This constraint corresponds to a single row in the rigidity matrix of the form 



b X a •••0 -- 



a X b 



Constraint (ii) requires a line to be parallel to a line in a fixed direction c 
translates into two rows: 



(ci,C2,C3). This condition 



V, 







(-C2,Ci,0) 




■ (C2,-Ci,0) 


...Q... 


(0, -C3,C2) 




1 (0,C3,-C2) 


...Q... 



The reverse direction of the proof of Theorem |16| requires us to assume that, except for the vectors v 
that are required to be zero in angular constraints, all other entries of the rigidity matrix are generically 
non-zero. Therefore, this representation of constraint (ii) is incompatible with the proof. 

However, since (ii) is expressing a line-line parallel constraint, we have an alternate description. Let a 
and b be two directions such that a x b = c. Then constraint (ii) can be expressed using two rows, each of 
which reduce to basic constraint (i): 



■0 



m 



c X a 
c x b 



a X c 


...Q... 


b X c 


...Q... 



For the basic blind constraints, let p be a point, p' its instantaneous velocity resulting from a twist and 
c ~ (ci,C2,C3) a direction vector. Then (Hi) expressed that p' must be orthogonal to c, and was associated 
to a single row 



0- (p : 1) V (c : 0) 



-(p:l)V(c:0) -0 



Constraint (iv) expressed that p' must be parallel to c and was associated to two rows 



...Q... 


(p:l)V(c2,-Ci,0,0) 


...Q... 


-(p: l)V(c2,-ci,0,0) 


...Q... 


...Q... 


(p:l)V(0,C3,-C2,0) 


...Q... 


-(p: 1)V(0,C3,-C2,0) 


...Q... 



As in constraint (ii), in constraint (iv) we are requiring that two vectors, p' and c be parallel. Again, 
there is an alternate description. Indeed, let a and b be two directions such that a x b = c. Then constraint 
(iv) can be expressed using two rows, each of which reduces to basic constraint (Hi): 



...Q... 


(p : 1) V (a : 0) 


...Q... 


-(p: l)V(a:0) 


...Q... 


...Q... 


(p : 1) V (b : 0) 


...Q... 


-(p: l)V(b:0) 


...Q... 
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We conclude that, generically, the entries in the three uj columns for basic angular constraints {i) 
and {ii) are non-zero, and the entries in the six v and uj columns for basic blind constraints {in) and 
(iu)are non-zero. Therefore, the body-and-cad rigidity matrix resulting from this alternative development 
will satisfy the definition of a (6, 3)-frame (as per Definition [3| . 

B Point-point coincidence constraints 

In this section, we provide some insight as to why we are unable to treat point-point coincidence constraints. 
If p = {p^ ,py ,p^) is the point of coincidence, then the constraint is infinitesimally expressed by 3 rows in 
the rigidity matrix of the form: 



...Q... 


(l,o,o,o,-p^p^) 


...Q... 


(-l,o,o,o,p^-p^') 


...Q... 


...Q... 


(o,l,o,p^o,-p-) 


...Q... 


(o,-l,o,-p^o,p-) 


...Q... 


...Q... 


(o,o,l,-p^p^o) 


...Q... 


(o,o,-l,p^-p^o) 


...Q... 



Since these rows do not conform to the pattern required by our proof technique, point-point coincidences 
are not addressed by the characterization. 

C Body-and-cad algebraic rigidity theory 

For completeness, we include the definitions from [s] required to express the geometry of a body-and- 
cad structure, from which the algebraic theory is derived. A cad graph (G, c) is used to represent the 
combinatorics of the constraints, where G = (V, E) is a multigraph with V — {1, . . . , n} and c : i? — > C is 
an edge coloring using colors C = {ci, . . . , C21} to denote the 21 possible cad constraints. This is related to 
the primitive cad graph presented in Section |2.1[ for each color q, we associate some number of primitive 
angular and primitive blind constraints as found in Table [l] (reproduced from [s]). The bi-colored primitive 
cad graph H = {V,E = BLiR) has a red edge for each primitive angular constraint and a black edge for 
each primitive blind constraint. 





point 


line 


plane 




angular 


blind 


angular 


blind 


angular 


blind 


point 


coincidence 





3 





2 





1 


distance 





1 





1 





1 


line 


coincidence 






2 


2 


1 


1 


distance 









1 


1 


1 


parallel 






2 





1 





perpendicular 






1 





2 





fixed angular 






1 





1 





plane 


coincidence 










2 


1 


distance 










2 


1 


parallel 










2 





perpendicular 










1 





fixed angular 










1 






Table 1: Association of body-and-cad (coincidence, angular, distance) constraints with the number of blind 
and angular primitive constraints. As an example of how to read the table, the last two columns (correspond- 
ing to plane) of row 3 (corresponding to coincidence under line) indicate that a line-plane coincidence 
constraint reduces to 1 angular and 1 blind primitive constraint. 
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The geometry of the constraint represented by each cad graph edge is given in terms of points, hnes, or 
planes affixed to the n rigid bodies; each hne or plane may be described by a point on a rigid body and 
a direction vector. Let Ei denote the set of edges of color a. Then we define a "length" function Li that 
associates to each edge in Ei the real- valued points, vectors and values (e.g., specified distance or angle) 
needed to describe the constraint; the full description of the family of functions Li is below. With these 
definitions, we can formally define a body- and- cad framework to be (G, c, L) = {G, c, Li, . . . , L21), the cad 
graph (G, c) together with the functions Lj. 

• Point-point coincidence: L\ : Ei ^ maps an edge e = ij to a point p so that it is constrained 
to lie on bodies i and j simultaneously. 

• Point-point distance: L2 : E2 ^ ffi'"* x M'' x M maps an edge e = ij to a triple (pi^pj, a) so that the 
point Pi affixed to body i is constrained to lie a distance a from the point pj affixed to body j. 

• Point-line coincidence: L^ : E^ ^ M.^ x (M'' x M:") maps an edge e = ij to a pair {pi, {pj,d)) so that 
point Pi affixed to body i is constrained to lie on the line (pj, d) afhxcd to body j. 

• Point-line distance: L4 : E4 ^ M.^ x {R^ x M.^) x R maps an edge e = ij to a triple (pi, {pj, d), a) so 
that point pi affixed to body i is constrained to lie a distance a from the line {pj, d) affixed to body j. 

• Point-plane coincidence: L5 : £5 ^ M'^ x (M'' x M^) maps an edge e = ij to a pair (pi, {pj,d)) so 
that the point pt affixed to body i is constrained to lie in the plane {pj,d) affixed to body j. 

• Point-plane distance: : Ef^ M"* x (R"^ x R^) x M maps an edge e = ij to a triple {pi, {pj,d). a) 
so that the point pi affixed to body i is constrained to lie a distance a from the plane {pj , d) affixed to 
body j. 

• Line-line parallel: Lj : Ef ^ M.^ x M.^ x M.^ maps an edge e = ij to a triple {pi,pj,d) so that the 
lines {pi, d) and {pj,d) affixed to bodies i and j, respectively, are constrained to remain parallel to each 
other. 

• Line-line perpendicular: Lg, : Eg, ^ {R"^ x R"^) x (M'^ x R^^) maps an edge e = ij to a pair 
{{pi,di), {pj,dj)) so that the lines {pi,di) and {pj,dj) affixed to bodies i and j, respectively, are con- 
strained to remain perpendicular to each other. 

• Line-line fixed angular: Lg : Eg ^ (M^ x M^) x (K^ x M^) x M maps an edge e = ij to a triple 

{{pi,di),{f>j,dj),a) so that the lines {pi,di) and {pj,dj) affixed to bodies i and j, respectively, are 
constrained to maintain the angle a between them. 

• Line-line coincidence: Liq : Eiq R^ x R^ maps an edge e = ij to a pair (p, d) so that the line 

(p, d) is constrained to be affixed to bodies i and j simultaneously. 

• Line-line distance: Ln : En — )• {R^ x R^) x {R^ x R^) x R maps an edge e = ij to a triple 
{{Pi,di),{pj,dj),a) so that the lines {pi,di) and {pj,dj) affixed to bodies i and j, respectively, are 
constrained to lie a distance a from each other. 

• Line-plane parallel: L12 : E12 — )■ {R^ x R^) x {R^ x R^) maps an edge e = ij to a pair ((pj, rf,), {pj, dj)) 
so that the line {pi,di) and plane {pj,dj) affixed to bodies i and j, respectively, are constrained to 
remain parallel to each other. 

• Line-plane perpendicular: L13 : E13 — > {R^ x M"^) x {R^ x R^) maps an edge e = ij to a pair 
{{pi,di), {pj,dj)) so that the line {pi,di) and plane {pj,dj) affixed to bodies i and j, respectively, are 
constrained to remain perpendicular to each other. 

• Line-plane fixed angular: L14 : £'14 — >• (R^ x R^) x (R^ x R^) x R maps an edge e = ij to a triple 

{{pi,di), {pj,dj),a) so that the line {pi,di) and plane {pj,dj) affixed to bodies i and j, respectively, are 
constrained to maintain the angle a between them. 
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Line-plane coincidence: L15 : £'15 (M"^ x M^) x (M^ x M^) maps an edge e = ij to a pair 
{{Pi,di),{pj,dj)) so that the hne {pi,di) affixed to body i is constrained to he in the plane {pj,dj) 
affixed to body j. 

Line-plane distance: Liq : Eiq — >• (M^ x R^) x {R^ x M^) x M maps an edge e = to a triple 
{{Pi^di), {pj,dj),a) so that the line (pi, di) affixed to body i is constrained to lie a distance a from the 
plane {pj, dj) affixed to body j. 

Plane-plane parallel: Ln : En — )• x x maps an edge e = ij to a triple {pi,Pj,d) so that 
the planes {pi,d) and {pj,d) affixed to bodies i and j, respectively, are constrained to remain parallel 
to each other. 

Plane-plane perpendicular: Lis '■ £-18 ^ {^'^ x x x maps an edge e = ij to a pair 
{(j)i,di),{pj,dj)) so that the planes {pi,di) and {pj,dj) affixed to bodies i and j, respectively, are 
constrained to remain perpendicular to each other. 

Plane-plane fixed angular: Lig : E'lg (M^ x M^) x (M^ x R^) x R maps an edge e = ij to a triple 

{{pi,di),{pj,dj),a) so that the planes (pi.di) and {pj,dj) affixed to bodies i and respectively, are 
constrained to maintain the angle a between them. 

Plane-plane coincidence: L20 : x maps an edge e = ij to the pair [p, d) so that the 

plane (p, d) is constrained to be affixed to both bodies i and j simultaneously. 

Plane-plane distance: L21 : E21 M.^ xR^ xR^ xM. maps an edge e = ij to a quadruple {pi,pj,d, a) 
so that the planes (pi, d) and {pj,d) affixed to bodies i and j, respectively, are constrained to have the 
distance a between them. 
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